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A  joint  Navy/NASA  experimental  investigation  was  conducted 
in  October  1983  in  the  Langley  Research  Center  VSTOL  Wind 
Tunnel  to  test  the  hypothesis  that  helicopters  with  unusually 
hub-weighted  radial  load  distributions  should  experience  a  more 
severe  first-harmonic  inflow  velocity  field  during  transition 
than  ordinary  helicopters.  This  report  presents  an  approximate 
analysis  of  the  experimental  results.  The  hypothesis  is 
strongly  supported.  Compared  to  the  rotor  configuration  with 
the  most  tip-weighted  load  distribution,  the  configuration 
with  the  most  hub-weighted  distribution  appears  to  have  experi¬ 
enced  an  approximately  50  percent  greater  first -harmonic  inflow 
at  a  50  percent  greater  critical  flight  speed. 
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ADMINISTRATIVE  INFORMATION 

This  analysis  and  the  Navy  participation  in  the  experiments  on  which  it  is 
based  were  funded  by  the  Independent  Exploratory  Development  Program.  The 
experiments  were  performed  in  the  4-  by  7-Meter  Wind  Tunnel  at  the  Langley 
Research  Center  courtesy  of  the  National  Aeronautics  and  Space  Administration 
(NASA)  . 
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INTRODUCTION 

In  late  1979,  flight  tests  were  undertaken  by  Kaman  Aerospace  Corporation 
on  an  H-2  helicopter  equipped  with  an  experimental  circulation  control  rotor 
(XH-2/CCR).  It  was  found  to  be  much  more  difficult  to  trim  the  helicopter  in 
pitch  at  low  speeds  than  anticipated;  in  fact,  the  flight  tests  were  eventually 
abandoned  with  no  flights  at  speeds  exceeding  the  critical  transition  speed. 

Although  a  number  of  problems  were  found  which  could  qualitatively  explain 
the  pitch  trim  difficulty,  it  was  hypothesized  that,  because  the  radial  load 
distribution  on  a  circulation  control  rotor  is  much  more  hub-centered  than  that 
on  an  ordinary  rotor,  the  transition  inflow  field  at  critical  speed  might  be 
more  severe.  Although  analyses  of  the  flights  completed  did  not  indicate  unusual 
Inflow  characteristics,  the  possibility  remained  that  unusual  effects,  including 
a  higher-than-normal  critical  transition  speed,  might  have  shown  up  had  further 
flights  at  higher  speeds  been  possible. 
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In  1983,  an  agreement  was  reached  between  the  David  Taylor  Naval  Ship 
Research  and  Development  Center  (DTNSRDC)  and  the  National  Aeronautics  and  Space 
Administration  to  conduct  experiments  in  the  Langley  Research  Center  (LRC)  VSTOL 
Wind  Tunnel  to  seek  evidence  concerning  possible  effects  of  radial  load  distri¬ 
bution  on  the  transition  inflow  field.  These  experiments  were  completed  in 
October  1983  by  a  joint  LRC/DTNSRDC  team  using  an  existing  DTNSRDC  circulation 
control  rotor  model.  This  report  presents  an  approximate  analysis  of  the  results. 

EXPERIMENTS 

The  rotor  model  was  80-in.  in  diameter  with  four,  5-in. -chord  blades.  These 
blades  were  judged  small  enough  to  yield  reasonably  valid  transition  data  in  the 
LRC  VSTOL  Wind  Tunnel  (Figure  1) .  The  model  was  equipped  with  a  pneumatic  valve 
system  capable  of  providing  blade  duct  mean  pressures  up  to  10  psig  and 
azimuthal  1-per-rev  harmonic  variation  of  the  pressure  of  amplitude  up  to  about 
60  percent  of  the  mean.  The  valve  characteristic  was  such  that  a  2-per-rev 
harmonic  variation  one-third  the  amplitude  of  the  1-per-rev  component  and  90  deg 
out  of  phase  with  it  also  resulted.  The  model  was  equipped  with  blade  duct 
pressure  transducers,  blade  root  flapping  moment  strain  gages,  and  azimuth 
reference  signal  generators.  A  more  complete  description  of  the  model  is  given 
in  Reference  1. 

The  model  was  mounted  on  a  four -component  DTNSRDC  load  cell  balance  designed 
to  measure  thrust,  pitching  moment,  rolling  moment,  and  yawing  moment  during 
hover  tests.  All  tests  were  performed  at  a  tip  speed  of  500  ft/sec.  The  thrust 
reading  of  this  balance  was  also  to  be  used  for  the  forward  flight  tests. 
Unfortunately,  the  balance  system  malfunctioned  during  early  hover  tests.  Enough 
data,  however,  were  acquired  before  the  malfunction  to  confirm  that  the  model  was 
functioning  properly  and  to  establish  an  approximate  correlation  between  the  first- 
harmonic  flapping  moment  strain  gage  measurements  and  the  actual  root  (hub-center) 
flapping  moment.  All  analyses  of  the  transition  data  then  had  to  be  performed  in 
terms  of  root  flapping  moment. 

A  tabulation  of  the  experimental  results  for  the  transition  inflow  analyses 
is  presented  in  Table  1.  In  addition  to  the  tabulated  data,  a  number  of  other 
data  channels  were  recorded  including  the  azimuth  angle  of  the  maximum  opening 


of  the  pneumatic  valve  (which  was  of  the  eccentric-circular-cam  type)  and  the 
output  of  blade  duct  pressure  transducers  located  at  the  tip  and  a  mid-radius 
position,  in  addition  to  the  "blade  root"  (10  percent-radius)  transducer  output 
used  in  the  analysis.  Comparison  of  the  first-harmonic  phase  angles  of  these 
various  quantities  with  each  other  and  with  the  recorded  cam  angle  reveals 
random  inconsistencies  in  the  phase  angles  on  the  order  of  ±10  deg.  Since  LRC 
was  unable  to  provide  on-site  reduction  of  the  high-speed  data,  the  unusually 
large  error  band  was  not  detected  in  time  to  initiate  efforts  to  identify  and 
correct  the  cause. 

In  addition  to  the  random  inconsistencies,  there  are  also  regular  incon¬ 
sistencies — notably  that  the  recorded  phase  angle  (i.e.,  azimuthal  position 
of  first-harmonic  peak)  of  the  duct  pressure  at  the  blade  tip  lags  that  at  the 
root  by  about  3  deg  on  average;  whereas,  due  to  the  time  lag  involved  in  near- 
sonic  transport  of  a  pressure  signal  from  root  to  tip,  the  blade  root  first- 
harmonic  pressure  peak  should  lead  that  at  the  tip  by  about  23  deg.  The  results 
of  the  analyses  suggest  that  the  recorded  tip  pressure  phase  is  more  nearly 
correct. 

AERODYNAMICS  PERFORMANCE  CODE 

The  performance  code  used  in  the  analyses  is  an  undocumented  code  developed 
by  the  author  on  an  HP-9836  desk  computer  for  quick-look  investigations  of  the 
effects  of  various  design  parameters  on  the  performance  of  circulation  control 
and  X-Wing  rotors.  The  code  employs  a  curve-fit  approximation  to  the  experi¬ 
mentally  determined  two-dimensional  characteristics  of  the  circulation  control 
airfoils  used  on  this  model  rotor,  including  Reynolds  number  and  compressibility 

effects.  It  also  employs  a  "distributed  momentum"  variation  of  the  approximate 

2 

inflow  equations  presented  by  Blake  and  White.  The  Blake  and  White  formulation 
can  be  written  as: 


where 


A*  =  V(A+y*4  -  m*2)/2  (2) 

X*  =  \i*/ (A+u*^  +  u*A  (3) 

The  "distributed  momentum"  variation  merely  replaces  the  quantity  A+  by  the 
quantity  A*/(  local  disk  loading) / (average  disk  loading)  . 

In  the  present  investigation,  a  further  modification  was  introduced.  The 
quantity  A*  was  replaced  by  the  quantity  A*  so  that  the  first-harmonic  in¬ 
flow  term  could  be  easily  varied  iteratively  (by  varying  K  m)  to  determine 
what  value  of  this  term  was  needed  for  best  correlation  with  the  experimental 
results.  For  purposes  of  this  investigation,  the  "rigid  rotor"  option  of  the 
code  was  used;  that  is,  elastic  deflections  of  the  blade  were  neglected. 

It  was  intended  to  validate  the  aerodynamics  performance  code  (and,  if 
needed,  make  empirical  adjustments  to  it)  by  comparing  calculated  to  measured 
thrust  and  hub  moments  in  hover  at  three  blade  angles  and  several  combinations 
of  blade  duct  mean  and  cyclic  pressure.  As  previously  mentioned,  only  one 
collective  angle  was  completed  before  the  thrust  balance  malfunctioned,  and 
only  a  few  cyclic  pressures  at  that  collective  angle  were  completed  before  the 
hub  moment  balance  began  to  behave  erratically. 

Correlation  of  the  few  valid  thrust  data  points  with  previous  measurements 
(unpublished)  and  with  the  calculated  behavior  is  shown  in  Figure  2.  It  is 
concluded  that  (1)  the  model  was  functioning  properly,  and  (2)  the  thrust  pre¬ 
dictions  from  the  aerodynamics  performance  code  are  satisfactory. 

The  few  valid  hub  moment  measurements  were  used  to  establish  the  approxi¬ 
mate  relationship: 
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Figure  2  -  Comparison  of  Measured  and  Calculated  Thrust 
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where  is  the  first -harmonic  amplitude  of  measured  flapping  moment  (measured 
by  a  strain  gage  bridge  at  the  12-percent  radius  station) .  From  this  it  is 
deduced  that : 


M,  =  1.25  Mr 
f  fe 

where  is  the  root  (hub  center)  flapping  moment.  This  relationship  was  employed 
in  all  subsequent  analyses . 

Correlation  of  the  measured  flapping  moments  with  calculated  behavior  is 
shown  in  Figure  3.  The  flapping  moment  predictions  from  the  aerodynamics  per¬ 
formance  code  are  also  satisfactory.  Consequently,  empirical  adjustments  to  the 
code  were  not  needed. 


ANALYTICAL  PROCEDURE 

No  inflow  quantities  were  directly  measured  in  the  experiment.  The  plan 
was  to  determine  for  each  experimental  data  point  the  value  of 


which,  when  substituted  for  Xj  in  Equation  (1),  would  yield  agreement  between 
the  measured  and  calculated  first-harmonic  root  flapping  moment  amplitude  and 
phase  (M^,!^).  The  calculation  would  have  employed  the  measured  values  of 
root  pressure  zeroth,  first-  and  second-harmonic  amplitude  and  first-  and  second- 
harmonic  phase. 

In  view  of  the  random  and  systematic  irregularities  in  the  pressure  phase 
measurements,  this  plan  was  modified  slightly: 

1.  It  was  assumed  that  the  best  first  estimate  of  actual  effective  root 
pressure  phase  was  a  value  leading  the  measured  phase  by  20  deg. 

2.  From  this  estimate,  random  variations  on  the  order  of  ±10  deg  can  be 
expected. 

3.  To  reduce  the  time  required  for  analysis,  it  was  decided  to  define 
"agreement"  by: 


BO  =  3.6 


O  EXPERIMENTAL 
□  CALCULATED 


o 


a.  Measured  Mf  cos  «  Calculated  Mf  cosi|»f. 

b.  Discrepancies  between  measured  and  calculated  Input  \p  should 
both  be  "reasonable"  In  light  of  the  above-discussed  data 
Irregularities. 

These  regrettably  soft  criteria  may  become  more  understandable  In  light  of 
the  example  presented  In  Figure  4.  Here  the  experimental  results  for  Run  61/Test 
Point  1  are  shown  with  the  calculated  results  for  various  combinations  of  input 
pressure  phase,  and  input  inflow  correlation  parameter,  Kgjtm.  Note  that  the 
"reasonable"  range  of  input  ^  (+1  deg  -  20  ±10  deg  -  -9  to  -29  deg)  brackets  the 
input  <J>p  of  -12  deg  required  to  produce  agreement  in  both  M^  cos  4*^  and  sin  ifi^. 
More  importantly,  note  that  the  value  of  required  to  produce  agreement  in 

M^  cos  alone  is  rather  insensitive  to  small  errors  in  either  ^  or  i|/  .  This 
insensitivity  to  errors  in  prevails  so  long  as  ip  is  near  zero  (or  180  deg)  , 
and  the  insensitivity  to  errors  in  \p^  prevails  so  long  as  ip^  is  near  180  deg 
(or  zero).  Fortunately,  a  goodly  fraction  of  the  data  were  taken  under  conditions 
meeting  both  of  these  criteria,  tending  to  strengthen  what  would  otherwise  be 
very  low  confidence  in  the  results  of  the  analyses. 

In  view  of  the  data  uncertainties  and  the  rather  time-consuming  calculation 
procedure  (about  10  min  per  calculated  point) ,  the  actual  analysis  procedure  did 
not  use  an  orderly  grid  of  calculated  points  such  as  illustrated  in  Figure  4. 
Rather,  from  a  first  guess  at  \p^  and  Kg^,  subsequent  iterative  estimates  (usually 
two  or  three)  were  made  until  the  result  was  judged  "close  enough,"  and  the  last 
estimate  was  recorded  in  Table  1.  Thus,  these  estimates  could  be  improved  upon 
by  further  analysis,  if  desired. 

The  second-harmonic  pressure  amplitudes  and  phases  are  not  recorded  in  Table 
1;  in  all  cases,  they  are  close  to  one-third  the  amplitude  of  the  first-harmonic 
and  90  deg  out  of  phase.  The  second -harmonic  component  is  included  in  the 
calculations.  The  effect  of  the  second-harmonic  component  is  to  make  the 
pressure  wave  broad/flat  on  top  and  narrow/peaky  on  the  bottom. 


RUN  61/1 
9  =  -3  deg 
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Figure  4  -  Effect  of  Variation  in  Input  Values  of  K£am  and  ipf  on 
Calculated  First-Harmonic  Components  of  Roof  Flapping  Moment 


RESULTS  AND  DISCUSSION 

The  results  are  summarized  as  plots  of  Aj  versus  p  in  Figure  5.  There  is 
a  systematic  variation  with  collective  blade  angle,  with  the  results  for 
0*0  deg  falling  in  reasonably  good  agreement  with  the  predictions  of  Reference 
1.  The  results  for  0  ■  -3  deg  and  0  =  -6  deg,  however,  fall  progressively 
further  above  the  predictions  of  Reference  1.  This  is  consistent  with  the 
original  hypothesis  that  the  first-harmonic  inflow  term  might  be  greater  for 
rotors  with  relatively  hub-weighted  radial  load  distributions. 

Figure  6  shows  the  calculated  blade  load  per  foot  at  the  advancing 
(i|»  *  90  deg)  and  retreating  *  270  deg)  blade  positions  for  two  thrust  co¬ 
efficients  at  each  of  the  blade  collective  angles.  At  0  *  0  deg,  blade  loading 
tends  to  increase  with  increasing  radius  over  most  of  the  blade  (as  it  does  for 
all  conventional  rotors);  at  0  =  -3  deg,  the  blade  loading  is  approximately 
constant  over  most  of  the  blade;  and  at  0  *  —6  deg,  blade  loading  decreases 
rapidly  with  increasing  radius  over  the  outer  half  radius. 

The  theory  of  Reference  1  is  a  rather  heuristic  one.  It  represents  the 

wing-like  component  of  rotor  lift,  as  determined  from  simple  momentum  theory, 

by  a  horseshoe  vortex  of  span  equal  to  rotor  diameter  and  estimates  Ax  from 

the  normal-to-disk  induced  velocities  of  this  vortex  system  on  the  centerline 

at  the  disk  leading  and  trailing  edges.  If  these  simple  assumptions  give  a 

correct  estimate  for  an  ordinary  rotor  (i.e.,  one  with  heavily  tip-weighted 

radial  load  distribution) ,  it  might  be  assumed  that  a  similar  estimate  using  a 

horseshoe  vortex  of  span  less  than  rotor  diameter  would  be  more  appropriate  for 

a  rotor  with  a  heavily  hub-weighted  radial  load  distribution.  Alternatively, 

this  speculation  can  be  expressed  in  the  form  that  such  a  rotor  should  behave 

(for  purposes  of  first-order  Inflow  estimates)  like  an  ordinary  rotor  at  the 

same  thrust  but  reduced  diameter;  or  like  an  ordinary  rotor  of  the  same  diameter 

but  greater  thrust  coefficient.  Such  an  estimate  is  achieved  by  plotting 
$  *  ±  * 

KA!  versus  Kp  ,  K>1,  where  Ai  (p  )  Is  taken  from  Equations  (1)  through  (3). 

Such  estimates  are  compared  to  the  results  of  Figure  5  in  Figure  7.  Values 
of  K  *  1.25  and  K  ■  1.53  are  found  to  correlate  the  0  ■  -3  deg  and  0  *  -6  deg 
data,  respectively.  (The  constants  1.25  and  1.53  were  selected  purely  to  fit 
the  data  and  not  from  any  theoretical  considerations.)  The  degree  of  data  fit 
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achieved  in  this  way  does  tend  to  support  these  heuristic  arguments  and  suggests 
that  some  further  simple  theoretical  explorations  seeking  to  connect  the  constant 
K  to  the  radial  load  distribution  in  a  quantitative  way  might  be  productive. 

Because  of  the  substantial  uncertainties  surrounding  the  experimental  data, 
it  is  necessary  to  question  whether  any  reasonable  systematic  error  could  produce 
the  trend  exhibited  in  Figure  7  artificially.  The  most  obvious  such  possibility, 
since  there  was  no  valid  measurement  of  thrust  during  the  transition  tests,  is 
the  possibility  of  a  systematic  error  in  the  calculated  thrust.  If  the  cal¬ 
culation  procedure  systematically  underestimated  the  thrust  at  0  =  -6  deg  and/or 
overestimated  the  thrust  at  0  =  0  deg,  an  artificial  difference  similar  to  that 
displayed  in  Figure  7  would  result.  However,  to  produce  the  magnitude  of 
difference  shown  in  Figure  7  would  require  that 


Calculated  thrust 


Actual  thrust  @  0  = 


@0=0  deg\  / 
=  0  deg  f  y 


Calculated  thrust  @  0  *  -6  deg' 


Actual  thrust  @  0  =  -6  deg 


be  equal  to  1.53  =  2.33.  This  is  regarded  as  extremely  unlikely. 

A  second  possibility  is  that  the  decision  to  base  the  analysis  on  root 
(hub-center)  flapping  moment  which,  in  turn,  was  estimated  as  1.25  x  (measured 
flapping  moment  at  12-percent  radius)  might  produce  a  systematic  error.  However, 
re-analysis  of  a  typical  data  point  at  each  0  was  performed  in  terms  of  the 


measured  12-percent  radius  flapping  moment  (with  no  additional  calibration  factor 

* 

applied).  The  result  was  to  reduce  the  estimated  value  of  by  about  3  percent 


in  each  case.  In  other  words,  there  was  no  effect  on  the  trend  displayed  in 
Figure  7. 

It  is  concluded  that,  although  significant  quantitative  errors  are  possible, 
the  trend  displayed  in  Figure  7  is  almost  surely  real,  and  rotors  with  unusually 
hub-centered  radial  load  distributions  do  indeed  experience  unusually  severe 


transition  inflow  fields. 
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Figure  5  -  Normalized  First-Harmonic  Inflow 
versus  Normalized  Advance  Ratio 
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Figure  5  (Continued) 
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Figure  7  -  Normalized  First-Harmonic  Inflow  versus  Normalized  Advance  Ratio 
Compared  with  Heuristically  Modified  Theory 
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